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We use analytical calculations and Monte Carlo simulations to determine the thermal fluctuation 
spectrum of a membrane patch of a few tens of nanometer in size, whose corners are located at 
a fixed distance d above a plane rigid surface. Our analysis shows that the surface influence on 
the bilayer fluctuations can be effectively described in terms of a uniform confining potential that 
grows quadratically with the height of the membrane h relative to the surface: V — (1/2)7/1^. The 
strength 7 of the harmonic confining potential vanishes when the corners of the membrane patch 
are placed directly on the surface (d — 0), and achieves its maximum value when d is of the order 
of a few nanometers. However, even at maximum strength the confinement effect is quite small and 
has noticeable impact only on the amplitude of the largest bending mode. 

PACS numbers: 



Fatty acids and other lipids are essential to every living organism. Because of their amphiphilic nature, they spon- 
taneously self-assemble into bilayer membranes that define the limits of cells and serve as permeability barrier to 
prevent proteins, ions and metabolites from leaking out of the cell and unwanted toxins leaking in In euokaryotic 
cells, membranes also surround the organelles allowing for organization of biological processes through compartmen- 
talization. In addition, biological membranes host numerous proteins that are crucial for the mechanical stability 
of the cell, and which carry out a variety of functions such as energy and signal transduction, communication, and 
cellular homeostasis 

An important aspect of biological membranes is that they are typically not free but rather confined by other 
surrounding membranes, adhere to other membranes, and attach to elastic networks like the cytoskeleton and the 
extracellular matrix. Several model systems with reduced compositional complexity have been designed to mimic 
biological membranes. These biomimetic systems include phospholipid bilayers deposited onto solid substrates (solid- 
supported membranes) Q, or on ultra-thin polymer supports (polymer-supported membranes) [J. With the aid 
of biochemical tools and generic engineering, supported membranes can be functionalized with various membrane- 
associated proteins _5j . Synthetic supported membranes with reconstituted proteins are increasingly used as controlled 
idealized models for studying key properties of cellular membranes [ff>] . They provide a natural environment for the 
immobilization of proteins under nondenaturating conditions and in well-defined orientation [7[ . Another attractive 
application of supported membranes is the design of phantom cells exhibiting well defined adhesive properties and re- 
ceptor densities Finally, biofunctional membranes supported by solid interfaces (semiconductors, metals, plastics) 
provide new classes of biosensors, diagnostic tools, and other biocompatible materials [1, Q. 

Theoretically, the thermal shape undulations of supported membranes have been addressed for various model 
systems. These model systems include: (i) membranes that adhere to surfaces under the action of a continuum local 
potential [Tol . [T]| . (ii) membranes pinned or tethered discretely to a surface [l^ [l3l [m , 15], and (iii) membranes 
supported by elastic networks of springs [H, [13, d, [H, [l^l • The investigation of the latter case is largely motivated 
by recent studies of simple cells (e.g., red blood cells and the lateral cortex of auditory outer hair cells), where the 
supporting cytoskeleton has a fairly well defined connectivity plj . One feature, missing in many previous theoretical 
studies, is the influence of steric ( "excluded volume" - EV) interactions between the membrane and its support on the 
elasticity and shape fluctuations of the membrane. When these interactions are considered (see, e.g., refs. [iTI. [2^). 
it is usually assumed that the disjoining potential due to the collisions between the membrane and the underlying 
surface decreases inverse-quadratically with the distance between them. This result has been originally obtained in 
a heuristic manner by Helfrich ^23J , and later was formulated more systematically by using a renormalization-group 
approach and computer simulations (see and refs. therein). An instructive way to understand this result is as 
follows: Consider a membrane that spans a planar square frame of area L^. The Helfrich energy (to quadratic order) 
for a nearly-fiat membrane in the Monge gauge is given by 



I. INTRODUCTION 




(1) 



where k the bending rigidity, and h the height of the membrane above the frame reference plane. Dividing the frame 
area into flo = {L/l)"^ grid cells of microscopic area P (where I is of the order of the thickness of the bilayer), and 
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introducing the Fourier transform of h (r) 



K = j2 f d'^rh{r}e"^'^ ; (f=^(ni,n2), ni,n2 = -^, ■•■ (2) 



the Helfrich energy takes the form 

^-j^Ef^l^^ (3) 
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from which (by invoking the equipartition theorem) one finds that 

(INR ^ ^g^. (4) 

where ks is Boltzmann factor, T the temperature, and q = \q\ ~ {2t: / L){n\+n\Y/'^ = {2i:/L)n. Using the last result, 
one readily finds that the mean-square fluctuation of the height increases quadratically with L 

q q 

Now, consider a membrane placed between two parallel walls positioned a distance d from each other, as shown 
schematically in Fig. [TJ The presence of the walls significantly suppresses the long wavelength thermal fluctuations 
of the confined membrane. At large length scales, we may assume that the net result of the confinement is that the 
membrane experiences an effective harmonic potential which can be introduced as an additional term in the Helfrich 
Hamiltonian: 



^ ^ (v2/i)' + 7/i2jdV, (6) 



where 7 is a constant which will be determined later, and h is measured from the mid-plane between the walls. In 
Fourier space the energy reads 

q 



and the spectrum of fluctuations is given by 



^\hq?) - J^^'l^l y (8) 



If 7 ^ k,{2tt/LY (which is always satisfied for sufficiently large L), then the mean-square fluctuation of the height is 
given by 



This result can be used for determining the value of 7. Each point of the membrane has equal probability to be found 
anywhere between the walls. Therefore, 



and by comparing Eqs. ([9]) and pO)) . we find that 



A' = ^, (10) 
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Due to the thermal fluctuations, the membrane collides with the confining walls and lose configurational entropy 
in these collisions. The walls will therefore experience a disjoining potential. Tracing over \hq\ in Eq.([7]) leads to the 
following expression for the Helmholtz free energy 



F 



27rL^kBT 



(12) 



where A is the thermal de-Broglie wavelength of a microscopic membrane patch of area P. The disjoining pressure 
between the walls is then calculated by using Eqs. ([5|)- P^ as follows: 



P 



1 dF _ 1 dFd-f _djA^ _ 6(fcsT)2 
From this result it follows that the effective disjoining potential per unit area [25j : 



V 



p{x)dx 



(13) 



(14) 



decays quadratically with d. 



II. MEMBRANE FLUCTUATIONS NEAR A SINGLE PLANE SURFACE 

Eq. fH]) for the disjoining potential has been derived for the case of a membrane fluctuating between two walls. 
Does this result also hold in the case of a membrane fluctuating near a single rigid wall? Consider a square membrane 
of linear size L with bending rigidity n whose four corners are held a distance d above a flat, impenetrable, surface 
(Fig. [2]). The height of the membrane relative to the underlying surface is denoted by the function h{x,y). In what 
follows, we shall assume that h{x, y) is periodic (with period L) along both x and y directions. The pressure due to 
the collisions between the fluctuating membrane and the surface must be repulsive. However, it is not a-priori obvious 
why p should be proportional to d~^, as predicted by Ea. ((T3|) . Moreover, it is not even intuitively clear whether 
this pressure should enhance or suppress the amplitude of the membrane thermal fluctuations. The pinning of the 
edges of the membrane and the EV interactions with the surface represent a combination of attractive and repulsive 
potentials whose net effect is not really well understood. A better understanding of this issue can be obtained by 
comparing the configurational phase space of our model system membrane with that of a freely fluctuating membrane. 
In the free membrane case, we consider the ensemble of configurations for which the spatial average of the height, 
h= (1/L^) J h{x, y) dxdy, is equal to some fixed arbitrary value c. Setting h — cis necessary to avoid multiple counting 
of physically equivalent configurations invariant under a vertical translation along the z direction. (Note the difference 
in notation between the spatial average A{h) = (1/i^) / dxdy A{h) which is calculated for a specific configuration, 
and {A{h)) which is the statistical mechanical average over the ensemble of all possible (distinct) configurations.) 
The phase space of free membrane configurations can be further divided into sub-spaces, where the height functions 
hi and /i2 of each two configurations included in the same sub-space can be related by hi(x,y) — h2{x + a,y + b) 
with < a,b < L. In other words, all the configurations in each sub-space can be transformed into each other by 
a horizontal translation in the x ~ y plane (see Fig. This transformation does not change h and, therefore, does 
not exclude (introduce) allowed (forbidden) configurations from (into) the phase space of configurations with h = c. 
Since the Helfrich energy [Eq.([T])] is invariant under translations, all the configurations within each sub-space have 
exactly the same statistical weight. The partition function which involves summing over all possible configurations 
can be presented as summation over the "sub-spaces" : 

Z= e "bt ^ ^ Q,e "bt ^ ^ e "bt , (15) 

configurations sub— spaces sub — spaces 

where is the number of configurations included in each sub-space. The last equality in the above equation can be 
understood as if each sub-space of the configurational phase space is represented by only one configuration with height 
function h{x, y) whose energy is given by by the sum of Helfrich elastic energy and an extra term that accounts for 
the "degeneracy" of the corresponding sub-space: 



-space 



= -kBT\nn + J |(v2/i)'dV. (16) 
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In the case of a free membrane, the number of configurations in each sub-space is obviously the same: Q. = Qq = (L/l)'^. 
(Note that for the purpose of counting the number of configurations, we henceforth assume that two configuration 
are distinct only if they are shifted by at least one grid cell of microscopic area P with respect to each other.) 

Let us repeat the above argument for our model system shown in Fig. [51 In this case, the mapping transformation 
between configurations belonging to the same sub-space involves two steps: (i) a horizontal translation in the x — y 
plane, and (ii) a vertical translation in the normal z direction which sets the height of the corners to be h{Q, 0) = 
h{0, L) = h{L, 0) = h{L, L) = d above the underlying surface (see Fig. [5] (A)). The Helfrich energy is invariant under 
these transformations. However, the vertical translation may lead to the intersection of the membrane with the surface 
and, therefore, to the exclusion of the configuration from the sub-space of allowed configurations (Fig. [4] (B)). The 
number of configurations left in each sub-space = il,{h{x,y),d) < Sip. Introducing the function G{h{x,y), d) < 1, 
we can write ft = J7oG(/i(a;, y), d) and rewrite Eg. p^ 

Wsub-space = -fci3rinr!o-fcBrinG(/i(x,2/),d)+ J ^(V^hfd^f. (17) 

Adding the term +kBT IhQq to Eq. (|17p allows us to replace the summation over "sub-spaces" back with summation 
over all the possible configurations of a "free" membrane (without a surface). The effective Hamiltonian of this "free" 
membrane is given by 

n = -kBT\nG{h{x,y),d) + J ^(V^/i)^rfV. (18) 

There is no EV term in this Hamiltonian (since the membrane is assumed to be free), but these interactions between 
the membrane and the surface are properly accounted for by the first term on the right hand side which quantifies 
the effect of the surface on the configurational entropy of the membrane. 

An interesting and unexpected result can be obtained in the d = limit, i.e., when the corners of the membrane are 
placed directly on the surface. In this limit, the transformation defined between configurations within each sub-space 
(see Fig. 131(A)) will almost always generate "forbidden" configurations that intersect the surface (Fig. 0] (B)). Only 
when the pinning point coincides with the global minimum of h{x,y), then the membrane will be positioned above 
the surface over the entire frame region. Therefore, each sub-space includes only one configuration: Q — 1 (neglecting 
the measure-zero set of configurations with multiple global minima), which means that G{h{x, y), 0) — fl/flQ = (l/L)'^ 
does not dependent on h{x, y). Substituting this result into Eq. pB)) . we find that a constant term was added to the 
Helfrich energy of the free membrane. Therefore, the statistical mechanical properties of the pinned membrane are 
identical to those of the free membrane and, in particular, its fiuctuation spectrum is also given by Eq.(|3]). 

How can we evaluate the function G{h{x,y),d) defined above? Let us consider again the mapping transformation 
between configurations belonging to the same configurational sub-space (Fig.[3](A)). This transformation changes the 
pinning point by translating the membrane both horizontally and vertically. The number of allowed configurations 
in the sub-space is determined by the number of points on the membrane which can be placed a height d above the 
surface without causing any part of the membrane to intersect the surface. The set of such possible pinning points 
includes all the points on the membrane for which h(x, y) — ft-min < where /imin is the global minimum of the height 
function. These points are located below the dashed horizontal line in Fig. [5l Denoting by Ap{h{x,y),d) < the 
total projected area associated with this set of possible pinning points, the function G{h{x, y), d) — Ap/L'^ represents 
the fraction of membrane points that satisfy the "pinning condition" h(x,y) — /imin < d. 

Let us introduce the height distribution function of the membrane, Ph(x.y){z)- For a given height function h{x,y), 
Ph{x,y){z)dz gives the fraction of membrane points for which z < h{x,y) < z + dz. The function G(h(x^y),d) is the 
cumulative distribution function associated with Ph{x,y){z)- 

G{h{x,y),d)^ Ph{x,y){z)dz. (19) 

We proceed by approximating p/j(j, j,)(z) by a Gaussian distribution function [26| 

1 ( {z-h?\ 
PHx.,y){z) - -=- exp —2 , (20) 

y/^T^^h{x,y) \^^h{x,y)/ 

where 

^^ix.,y)^ih~hy = (21) 
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Using Eq.dlQl) in Eq.(IT9l), we find : 

G{h{x,y),d)^l 



1 + erf 



Klin - h + d 



V2A 



h(x,y) 





1 




)] 




1 + erf 1^ 




^ 2 





V2A 



h{x,y) 



(22) 



where erf(x) = (21 ^pn) J^^ e^" du is tlie standard error function, and a = (h — Ka\-a) I \f^Ah(x.y) > 0. The function 
G{h{x,y),d) given by Ea. ((22|) satisfies the boundary condition that G{d oo) 1. The value of a can be set by 
imposing the other boundary condition (see discussion above) that G{d = 0) = (l/L)'^, which gives 



a 



erf 



2P^^ 
L2 



(23) 



Using the Fourier representation {hq} of the function h{x,y) and Eqs. (fTS]) and ((22)) . we find that the statistical 
mechanical properties of the pinned membrane can be derived by considering a free membrane whose thermal behavior 
is governed by the Hamiltonian: 



n = -fcsTln 



1 



rf 



V2A 



11 
L2 



(24) 



where a is given by Eq. ((23| . The dependence of the first term on the right hand side of Ea. ((24| on {hq} is contained 
in the variable A^i^^ y) (see Eq. ([2T|) ). The spectral intensity of the pinned membrane can be calculated ed by using 
the equipartition theorem: 



(25) 



Introducing the variable d* = —a 

. dn 



d/V2A 



h{x,y) 



the right hand side of Eq. ((25)) can be written as follows 



'dha 



dd* 



Lj nq'\hq\' - ksT ( 1 :^^hq 



) tiq^\hq\^ + ksT 



0F[l + erf (d*)] / \ 



h(x,y 



dAh(x^ 




'2 e 



A3 

^h(x,y] 



I V TT [1 + erf (d*)] 




(26) 



Eg. ipS)) represents a set of linear equations (one for each Fourier mode q 0). These equations are coupled to each 
other through the mean-square height fluctuation A^^ y-^ (see Eq. ipT]) ) appearing both explicitly on the second term 
in the curly brackets, as well as in the definition of the variable d* . For both d ^ oo (free membrane) and d — 
(membrane pinned directly to the surface) the second term in the curly brackets vanishes and, therefore, Eq.Q which 
describes the fluctuation spectrum of a free membrane is recovered in these two limits, as argued above. For finite 
values of d, a further approximation can be made by replacing the spatial average A/j^^, j,) with Aq, the ensemble 
average over free membrane configurations (see Eq.([5])). This approximation leads to the decoupling of the set of 
equations ([26)1 and yields the following result: 



kBTL"^ 



knT- 



W) \ V [l+crf(d*)] y 



(27) 



where within the approximation of replacing A/j^^, j,) with Aq, we also set d* ~ —a + d/\/2Ao. Eq. (l27|) has the same 
form as Eq. ([5]) which describes the power spectrum of a membrane fluctuating under the action of a uniform harmonic 
potential of strength 



7off 



keT 




(28) 
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From Eq. p7|) we conclude that the thermal height fluctuations of the pinned membraned are attenuated compared to 
the fluctuations of a free membrane. The mean square fluctuation amplitude of a mode with wavevector q = 2Tm/L 
is reduced by a factor 



In 



Kq 



(l^(7cff)P) 



(29) 



where 



/7cff\ / L_ 




1 + erf ( —a 



(30) 



is a dimensionless number that governs the crossover between the regimes of damped {n"^ n^) and free (n^ 3> n^) 
thermal fluctuations. Values of n'^ are plotted in Fig.[n]for different values of d and for k = lOksT and L = 101 ~ 50 nm 
(a = 1.645 - see Eq.^). The maximum value of - 0.34 is obtained for d - 0.025 - 0.04L - 1.25 - 2 nm. This 
value of ^ 0.34 is too small to have any noticeable effect on the spectrum of thermal fluctuations, except for 
the largest mode (n = 1) whose square amplitude is suppressed by a factor of Ii ~ 1/1.34 = 0.75. In comparison, 
the square amplitudes of second and third largest modes are reduced by only factors of ^ 0.92 and I2 ^ 0.98, 
respectively. We also observe from Fig. [S] that for smaller values of k, the maximum values of occurs at larger 
values of d/L. This behavior is anticipated since the smaller k the larger the amplitude of the thermal fluctuations 
and, therefore, the greater the range of steric repulsion between the membrane and the surface. However, even for 
very soft membranes with k = SksT, the maximum is still reached at d/L < 0.1, i.e., only a few nanometers above 
the surface. 



III. COMPUTER SIMULATIONS 



One of the important recent advances in soft-matter simulations is the development of coarse-grained (CG) bilayer 
membrane models in which the membranes are simulated without direct representation of the embedding solvent 
[23, HI]. These implicit solvent ("solvent free") CG models require modest CPU and memory resources and, therefore, 
can be used for simulations of mesoscopically large membranes over long enough timescales to address experimental 
reality. Here, we use an implicit solvent CG model to test the validity and accuracy of the analytical predictions 
discussed above. A snapshot from the simulations is shown in Fig. [71 Each lipid molecule is represented by a short 
string of three spherical beads of diameter a, where one of the beads (depicted as a dark gray sphere in Fig. [7]) 
represents the hydrophilic head group and two beads (light gray spheres in Fig. [7]) represent the hydrophobic tail of 
the lipid. The details of the model and of the molecular simulations are given in refs. [2^ and including the 
description of a new Monte Carlo scheme (Mode Excitation Monte Carlo) which has been applied to accelerate the 
relaxation of the five largest Fourier modes corresponding to n = 1, V2, 2, -\/5, -\/8. The simulated system consists 
of a bilayer membrane of 2000 lipids and a surface (whose boundary is indicated by a thick black line) located at 
z = below which the beads cannot be found. The center of the head bead of one of the lipids in the lower leaflet 
(indicated by a black sphere appearing at the front of FigO is held fixed at r = (x, y, z) = (0, 0,d + a/2). Note that 
in our simulations, d is defined as the distance between the surface and the bottom of the particle whose position is 
held fixed. The simulations were carried out on the "high performance on demand computing cluster" at Ben Gurion 
University. For each value of d, the simulation results appearing below are based on 16 independent runs, each of 
1.2 X 10^ Monte Carlo time units. The first 2 x 10^ time units of each run were discarded from the statistical analysis. 

The simulations were conducted in the constant surface tension ensemble, at vanishing surface tension. We set the 
bead diameter, ct ^ 5/6 nm, so that the bilayer membrane thickness is Z ~ 6(T ~ 5 nm, and the (mean) linear size 
of the bilayer L = (36.4 ± 0.1) u ~ 30 nm (for all values of d). The spectrum of the membrane height fluctuations 
is plotted in Fig. [5] for both d = and d = 00 (i.e., for a freely fluctuating membrane). The computational 
results fully confirm our analytical prediction that the power spectra in these two cases are identical. By fitting the 
computational results to the asymptotic form (for small values of n) Eq. (|4]), we obtain that the bending rigidity of 
the membrane k ~ (7.8 ± 0.2) ksT. 

For d ^ we expect the power spectrum of the membrane to be "almost" identical to the power spectrum of the 
free membrane. More precisely, Eq. predicts that the mean square fluctuation amplitudes of the modes will be 
reduced by a factor J„ which, except for the longest n — 1 mode, is very close to unity. This prediction is very well 
supported by our computational results which are summarized in Table [H The table gives the values Ii and 
(corresponding, respectively, to the largest and second largest Fourier modes) for different values of d. Within the 
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statistical accuracy of our simulation results, we found no change in the fluctuation amplitudes of all the other modes 
corresponding to wavenumbers n > 2. In order to evaluate the quantitative agreement between the computational 
data presented in Table [J and the above analytical theory, we use Eqs. ([5]), (|23| . (|29| . and ([30]) to calculate the factor 
Ii for different values of d and for the set of parameters relevant to our simulations: k — 7.8kBT, L = 61 = 36(T. The 
results of the calculation along with our computational results (Table [l| are plotted in Fig.[9l As can clearly be seen 
in the figure, the agreement between the analytical and computational results is quite good. This agreement lends 
support for the validity and accuracy of our theoretical analysis of the fluctuation spectrum. 



IV. THE DISJOINING PRESSURE 

The fact that surface influence on the bilayer fluctuations can be effectively described in terms of a uniform harmonic 
potential, does not imply that the disjoining pressure between the surface and the membrane follows Eq. ()13|) . The 
pressure can be related to Hamiltonian (j24|l through the following equation: 



j_/dn\ 

' L^\dd / 



(31) 



from which we readily derive that 



L2 



V 



-a+ 



V2A, 



1 + erf ( -a 




(32) 



Following the approach described in section |TT1 the thermal average in Eq. (|32p can be approximately evaluated by 
replacing A^(^ with Aq (Eq.Q), which gives the following expression: 




1 + erf ( —a 



= P 



V2A0, 



[1 + erf {d* 



pH{d*), 



(33) 



where p* = {ksT nfl"^ j , d* = -a + d/V^Ao, and C = {2tt)^{2/6.03tt)^/^ = 12.83. The scaling function H{d*), 
which is plotted in Fig. [TOl decreases monotonically with d* . At large distances {d — > 00) the pressure decreases as 
~ exp[— (d/L)^]. The maximum pressure occurs when the membrane is in direct contact with the surface (d — 0). 
The contact pressure, however, does not diverge but rather reaches the following finite value p{d = 0) = p*H(—a) ~ 
p*^/TTCa. Using Eqs. (|23p and p3p . one can easily realize the contact pressure diminishes with the size of the 
membrane patch: 



p{d ^0)= p*^C ■ erf"^ [l - (2/7^^)] 



VkBTK ( L 



(34) 



V. CONCLUDING REMARKS 



In this paper, we analyzed the influence of EV volume effects on the statistical mechanical properties of supported 
membranes. Using analytical calculations and Monte Carlo simulations, we determined the fluctuation spectrum of 
a bilayer membrane patch of a few tens of nanometer in size whose corners are located at a fixed distance d above a 
plane rigid surface. We found that the surface has influence on the fluctuation spectrum of the membrane only when 
the pinning distance d is sufficiently small {d < O.IL - see Figs. [S] and [S]). This result can be easily understood given 
the fact that the amplitude of the height fluctuations of a free membrane patch satisfies Aq ~ 0.02L (see Eq. ([S])) 
and, obviously, the membrane hardly collides with the underlying surface when d ^ Aq. At small distances {d ~ Aq), 
the surface influence on the bilayer fluctuation spectrum resembles that of a uniform harmonic conflning potential 
of the form: V — {l/2)^cs{d)h'^ ■ Both analytically and computationally we find that the strength 7off(d) of the 
effective harmonic potential is extremely small and has noticeable impact only on the amplitudes of the very largest 
fluctuation modes. More remarkably, the confinement effect vanishes (7cff(d) = 0), when the membrane is brought 
into direct contact with the surface (d = 0). This unexpected and counterintuitive result can be explained by the fact 
that the thermal motion of the membrane is not really confined within a finite spatial region (as in a stack of bilayer 
membranes) but only restricted on one side. Therefore, the primary effect of the collisions with the underlying surface 
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is to push the membrane "upward" rather than to suppress the ampHtude of the fluctuations. This is also the reason 
why the disjoining pressure does not diverge when d vanishes. 

The significance of our findings should be considered in light of previous theoretical attempts to quantify the steric 
effects between membranes and underlying supported interfaces. Some of these studies describe the wall-membrane 
pressure by means of Ea. (|13|) . which we have shown to be irrelevant for this problem. One particular problem that 
should be reanalyzed in light of our new results is the theoretical interpretation of the fluctuation spectra of red 
blood cells [sT'l. The plasma membrane of red blood cells is attached to a triangulated network of flexible spectrin 
proteins with mesh size ^ ~ 60 — 100 nm. The spectrum of red blood cells fluctuations was analyzed in terms of the 
Helfrich Hamiltonian with both curvature and scale dependent surface tension terms, where the latter term originates 
from the coupling to the cytoskeleton. Using a Gaussian network model, Fournier et al. [13, [11] showed that the 
effective surface tension exhibits a steep crossover from a vanishingly small value at length scales smaller than ^ to 
some finite value at scales larger than ^. Gov et al. [l6l [2l| argued that, in addition, a uniform harmonic potential 
term must be introduced in the Helfrich effective surface Hamiltonian, which accounts for the confinement effect due 
to the steric repulsion between the spectrin and the bilayer. Our statistical mechanical analysis partially supports 
this phenomenological argument. On the one hand, our Eq. ()27p can be interpreted in terms of a uniform harmonic 
potential that acts on the membrane. On the other hand, our estimation of the strength of the effective harmonic 
potential makes it questionable whether the origin of it can be attributed to EV interactions (between the bilayer 
and the spectrin) alone. It seems more likely to relate this additional confinement term to the junctional complexes 
(of short actin filaments, globular band 4.1, and other proteins) which connect the membrane to the cytoskeleton 
and restrict the membrane height fluctuations around the points of attachment. We thus speculate that, just like the 
surface tension, the strength of the effective harmonic potential 7eff must also be scale dependent. At length scales 
below the mesh size, we expect the value of 7cff to be governed by EV effects and, therefore, to be extremely small. 
Above the mesh size, the strength of the harmonic confinement will be determined by the strength of the periodic 
pinning of the membrane to the cytoskeleton which, presumably, result in a larger value of 7efi. It should be stressed 
here that there is currently no proof (or even a reasoned argument) that the long wavelength fluctuations of red cells 
are indeed harmonically confined. One should also bear in mind that on the scale of the mesh size of the network, 
the problem is quite intricate and issues such as connectivity defects and the motion of the protein anchors must be 
properly addressed. 

Discussions with Nir Gov, Thorsten Auth, Sam Safran and Phil Pincus are gratefully acknowledged. 
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d 


h 




0.5(7 


0.89 (5) 


0.99 (4) 


a 


0.87 (6) 


0.97 (4) 


1.5a 


0.84 (4) 


0.94 (4) 


2a 


0.85 (4) 


0.95 (4) 


2.5a 


0.90 (6) 


0.96 (4) 


3.5a 


0.94 (4) 


0.99 (4) 


4.5a 


1.00 (5) 


1.00 (4) 



TABLE I: The factors 7i and (see Eq.([29}) by which the mean square fluctuation amplitudes of the largest (n = 1) and 
second largest in — \f2) modes are reduced as compared to the square fluctuation amplitudes of a free membrane. The height 
d denotes the distance between the bottom of the flxed spherical bead and the underlying surface. 



d 




I 

^ 

FIG. 1: A fluctuating membrane conflned between two walls which are separated from each other by a distance d. The height 
of the fluctuating membrane, h, is measured from the mid plane between the walls (— d/2 < h < +d/2). 




FIG. 2: A square membrane fluctuating above a flat, impenetrable, surface. The function h{x,y) denotes the height of the 
membrane above the underlying surface. At the four corners of the surface h{0, 0) = h{0, L) = h{L, 0) — h{L, L) = d. Outside 
the frame region, h{x, y) is deflned by periodic extension. 
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hi h2 




L 

FIG. 3: Configurations that arc represented by the height functions hi{x,y) (sohd curve) and h2{x,y) (dashed curve) belong 
to the same sub-space of configurations if these configurations are invariant under translation in the x-y plane, i.e., hi{x,y) = 
h2{x + a,y + b) with < a,b < L. 




FIG. 4: (A) The transformation between the configurations represented by the height functions hi{x,y) and h2{x,y) involves 
(i) a horizontal translation in the x — y plane and (ii) a vertical translation in the normal z direction that sets the corner at 
hi{0, 0) — /i2(0, 0) — d. (B) Configurations such as the one represented by the function hs intersect the underlying surface and, 
therefore, should be excluded from the phase space. 
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FIG. 5: The membrane can be pinned to the surface without intersecting it only at points for which h{x,y) — hmin < d, i.e., 
the membrane points located below the horizontal dashed line in the figure. Specifically, for d = the only possible pinning 
point is at the global minimum of the function h. 




FIG. 6: The dimensionless parameter (see definition in text - Ea l30|l as a function of the height of the pinning points from 
the surface. The values have been calculated for a membrane of linear size L = lOZ with bending rigidity k — lOfcaT (solid 
line) and k = 3KbT (dashed line). 
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FIG. 7: Equilibrium configuration of a membrane consisting of 2000 lipids. Each lipid is represented by a trimcr of one 
"hydrophilic" bead (dark gray sphere) and two "hydrophobic" beads (light gray spheres) . The membrane is fluctuating above a 
plane surface (frame indicated by a thick black line) , while the position of the center of one of the hydrophilic beads (appearing 
at the front of the figure and indicated by the black sphere and an arrow) is fixed at r = (x, y, z) = (0, 0, d + cr/2). 




2 



FIG. 8: The fluctuation spectrum of a membrane of N = 2000 lipids. The results for d = (membrane pinned directly to the 
surface) are shown by small solid circles. These results are essentially identical to those obtained from simulations of a free 
membrane (d — oo), which are represented by larger open circles connected with a dashed line. The solid line indicates the 
asymptotic (Iftql^) ~ power law. 
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